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In this paper, we analyze the clustering of galaxies using a modified Newtonian potential. This 
modification of the Newtonian potential occurs due to the existence of extra dimensions in brane 
world models. We will analyze a system of galaxies interacting with each other through this modified 
Newtonian potential. The partition function for this system of galaxies will be calculated, and this 
partition function will be used to calculate the free energy of this system of galaxies. The entropy and 
the chemical potential for this system will also be calculated. We will derive explicit expression for 
the clustering parameter for this system. This parameter will determine the behavior of this system, 
and we will be able to express various thermodynamic quantities using this clustering parameter. 
Thus, we will be able to explicitly analyze the effect that modifying the Newtonian potential can 
have on the clustering of galaxies. We also analyze the effect of extra dimensions on the two-point 
functions between galaxies. 


I. INTRODUCTION 


Motivated by the existence of extra dimensions in string theory, brane world models propose that our universe is 
a brane in a higher dimensional bulk 0,i- The extra dimensions can be either compactified or infinite, and the 
number of such extra dimensions along with the field content of the theory depends on the details of the model that is 
being investigated. However, common feature of all of these different brane models is that the standard model fields 
are confined to the four dimensional brane. The gravitons are not confined to the brane and can propagate into the 
higher dimensional bulk. This explains the weakness of the gravitational field [I|, 0 BjS- In fact, it is hoped that 
these models can also provide a solution for the solution of the hierarchy problem [ 2 , Si 1 Jt 3 ■ may be noted that 
the Newtonian potential gets modified at large distances due to the propagation of gravity into the bulk Q. These 
corrections are parametrized by Yukawa potential Q. Thus, because of the existence of extra dimensions, the original 
Newtonian potential is modified to Q 

$ = $jv(l+ae^“), (1) 


where a, A depend on the details of model, and Vm is the original Newtonian potential, 

Gm^ 
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However, in the Randall-Sundrum model the modification to the Newtonian potential is parameterized by a power 
law deformation of the usual Newtonian potential 0,i 


$ = (^1 + 

The action for the Randall-Sundrum model is given by Q 
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Here R 5 is the five dimensional scalar curvature, 35 and are determinants of the five dimensional and the four 
dimensional metrics, respectively. The M denotes the a suitable constant formed from the five dimensional Planck 
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mass. The brane tension is denoted by Vbrane- The four dimensional metric is induced on the hyper-surface y = 0. 
In the weak field approximation, this model gives a power law deformation of the usual Newtonian potential (3j. 

Unlike many other brane world models, the Randall-Sundrum model allows for the existence of infinite extra 
dimensions. It is this power law deformation of the usual Newtonian potential which we will use for our analysis in 
this paper. It is worth mentioning that there are different app roaches where Newton’s law of gravity get modified. 
One of these approaches is the non-commutative geometry [iflllli- Besides, some modifications of Newton’s law of 
gravity have been found due to minimal length in quantum gravity [T^ . In addition, a modification of Newton’s law 
of gravity has been predicted in entropic grav ity |13j . Besides, another sources of modifications to Newton’s law of 
gravity have been investi gat es in f(R) theories |14| . It may be noted that even a modification based on entropic force 
produces similar results |15l Il6l| . Even though it is related some minimal length deformation, the important thing 
to note here is that this modification of Newtonian potential can be motivated from various physical phenomena. 
So, in this paper, we will use the modified Newtonian potential based on Randall-Sundrum models for analyzing the 
clustering of galaxies. 

The modification of the Newtonian potential can occur only at very small or very large distances, at which the 
Newtonian theory of gravity is not tested. In fact, it has been observed that the rotation curves of galaxies do not 
fit the theoretically predicted curves, and this has motivated the development of theories like the MOND [i3[i3 and 
MOG [Il,[2^. So, we expect the deviation from the Newtonian potential to become effective at galactic scale. In this 
paper, we will study a system of of galaxies interacting gravitationally with each other. This system is large enough 
the modihcation to the Newtonian potential to become significant. We will analyse the effect this modification of the 
Newtonian potential has on the clustering of galaxies. In our analysis, we will effectively treat each galaxy as a point 
particle. However, we will incorporate an random parameter, which will take into account the fact that galaxies do 
actually have an intended structure. As we will be dealing with a large number of galaxies, and individual galaxies 
will be effectively treated as point particles, we can apply the methods of statistical mechanics and calculate the 
partition function for this system l25l - l^ . The strength of interaction will depend effectively on the average 

density of the universe [2lj. 

A critical physical assumption in this analysis is that the local dynamical timescale in an over dense region varies 
faster than global gravitational time scale. This assumption makes it possible for gravitational clustering to evolve 
through a sequence of quasi-e quil ibrium state [l^. In fact, using this assumption the cosmological many body partition 
function has been calculated |26l. . This was done by using using an ensemble of co-moving cells, and each of 

these co-moving cells contained galaxies which interacted with each other. It was demonstrated that if the size of the 
cells is smaller than the particle correlation length, then each member of this ensemble is correlated gravitationally 
with other cells. Furthermore, the correlations within each cells is greater than correlations among different cells, so 
that extensive is justified as a good approximation [^. l30l|. In this paper, we will use this formalism for analyzing a 
system of galaxies interacting through a modified Newtonian potential. 


II. GRAVITATIONAL PARTITION FUNCTION 

In this section, we will calculate the partition function for a system of galaxies interacting though a modified 
Newtonian potential. We consider a large system, which consists of an ensemble of cells, all of the same volume V, 
or radius i?i and average density p. Both the number of galaxies and their total energy will vary among these cells, 
thus these are represented by a grand canonical ensemble. Now the galaxies within the system have gravitational 
pairwise interaction generated by the modified Newtonian potential. It will be further assumed that the distribution 
is statistically homogeneous over large regions. 

The general partition function of a system of N particles of mass m interacting through the modihed gravitational 
potential with potential energy is $, can be written as 


= ^^ 3 ^ J X exp(^- ^ ^ + $(ri,r2,r3,... ,rjv) T 


(7) 


where pt is the momentum for different galaxies and T is the average temperature. It may be noted that a sim¬ 
ilar analysis has been done for galaxies interacting though the usual Newtonian potential [^. Here N\ takes the 
distinguish-ability of classical particles into account, and A is the normalization factor which results from integration 
over momentum space. Now integrating the momentum space, we obtain the following expression 
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where Qn{T,V) is the configuration integral 


Qn{T, t") 
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The gravitational potential energy $(ri,r2 ,... j'tat) is a function of the relative position vector Vij = — rj\ and 

is the sum of the potential energies of all pairs. Even for a system of galaxies in brane world models, the potential 
energy $(ri, r2 ,..., r^) can be expressed as 

$(ri,r2,...,riv) = JZ 

In fact, we can use a two particle function — 1, such that it vanishes in absence of interactions, and is 

non-zero only for interacting galaxies. 

It has been demonstrated that for galaxies interacting through a usual Newtonian potential, that the configurational 
integral can be expressed in terms of the function [2^. Repeating this analysis for galaxies interacting through 
modified Newtonian potential, we obtain 


Qn{t, V) 
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The interaction potential energy between galaxies in brane world model with large extra dimensions, can be written 
as 0,0 


4>i,j — 



( 12 ) 


where r » and is considered as a typical length scale at which the correction dne to the infinite extra 

dimensions become dominant. It may be noted that for point masses the partition function for galaxies interacting 
through usual Newtonian potential diverges at = 0. This divergence occurs due to the assumption that the 
galaxies are point like objects. However, this divergence can be removed by taking the extended nature of g alaxies 
into account by introducing a softening parameter which takes care of the finite size of each galaxy (2^, [2^, [23. Thus 
by incorporating the softening parameter the interaction potential energy between galaxies in brane world models can 
be represented as 
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where e the softening parameter which is related to the extended nature of galaxies. Now we can write 
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Expanding the function fij, we obtain the following expression. 
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Now we obtain (5i(T, V) = V hy evaluating the configuration integrals over a spherical volume of radius Ri. Similarly, 
evaluating the configuration integrals for Q 2 (T, V ), we obtain 


Q 2 {t,v) 




1 ,Gto2 /■«! 


n=0 


nl 


T 


r 

Jo 


(7’2 _j_ ^2'jnj2 



ki 

(4 + e2)V2 


n 

dr. 


(16) 






























4 


For large values of r, we can neglect the higher terms corresponding to _|_g 2 )i /2 ? a-nd thus, we can write 
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g.(T,h^) = 4.i/f;l(^rx 
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Evaluating this integrals, we obtain 


Q2(r,E) = E2(i + ^(a„x"+/3„x")), 


n —1 


where x = 3/2/3pT The values of a„ and j3n are given by 
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where E^|, is the hyper-geometric function. Now as i?i ^ p we can write 


3Gm^ 3Gm? 3 




2 RiT 2p-i/3T 2 

Using scale invariance, p —>■ X~^p, T —)> A“^T and r —)> Ar, we obtain 


-{Gm^fpT-^ = PpT-^ 
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where /3 = ^{Grri^)^. Thus, we can write 

00 

Q2(T,U) = U2(l + ^(a„ + /3„)^”), 

n—1 

where x = (3pT~^. By similar procedure, we obtain configuration integral an expression for general term as 

00 

Qn{T, U) = U'^(1 + ^(a„ + Pn)x^f~\ 

n—1 

Thus, we can write 

00 

Q^{T,V) = U3(l + ^(a„ + /3„)x")^ 

n—1 
00 

Qa{T,V) = V\\FY,{ar, + Pn)x^f. 

n—1 

Hence, the gravitational partition function for a system of galaxies in the brane world model, can be written as 

Zn{T,V) = ±(“)^^/V^(l + + (26) 


( 22 ) 

(23) 

(24) 

(25) 


It may be noted that the weakly interacting system can be represented by n = 1, and the strongly interacting systems 
can be represented by taking higher values of n. 
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III. THERMODYNAMICS OF GALAXIES 


Now we will use the partition function for the system of galaxies interacting in a brane world model to calculate 
relevant thermodynamic quantities. We can define for different values of n as 


Br, = 


Y^^=in{an + Pn)x‘^ 


(27) 


This is the general clustering parameter for gravitationally interacting system of galaxies in brane world models. We 
can express various thermodynamic quantities in terms of this general clustering parameter. Thus, we can write the 
Helmholtz free energy F = —TlnZN{T,V) as 


n—1 ' 


(28) 


The entropy S can now be calculated from this Helmholtz free energy, 
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where S'o = §77 + §77 In . The internal energy U = F + TS and of a system of galaxies can now expressed as 


17 = -Ar(l- 2 B„). 

Similarly, we can write the pressure P and chemical potential ^ can be expressed as 
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The probability of finding N particles can be written as 
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and z is the activity. Thus, for a system of gravitationally interacting system of galaxies in brane world models, we 
can write 
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Now the grand partition function for this system of gravitationally interacting system of galaxies in brane world 
models,can be written as 


lnZG=^=iV(l-i3„), 


(37) 


and the distribution function can be expressed as 
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This is the general distribution function for any order of terms. Thus, we can write for n = 1, which is the case of 
weak interactions, and n = 2 and n = 3, 
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Similarly, we can calculate the higher order distribution functions by taking higher order terms. More the number of 
terms more clustering in the system. Thus, by including more terms we may have the distribution function of strongly 
interacting system of galaxies in the brane world models. 


IV. COMPARATIVE STUDY 


We will be comparing our results with our previous work. As per the previous work the clustering parameter with 
first term is expressed as [1^ 


b = 


aix 

1 + aia; 


(40) 


Now after inclusion of brane corrections the clustering parameter can be expressed as 
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The relation between B and b is obtained and we express it as 
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where for n = 1 


fczi^(3/2,(; + l)/2;5/2;-i?7e2) 
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where we have put 
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This ratio can be calculated for different values of 1. e.g for Z = 1 a = 0.29, I = 2 a = 0.103, I = 3 a = 0.046, I = 4 
a = 0.025 assuming ki = e, 2e We have made a graphical comparison of distribution function F(N) for different values 
of I to see the effect of first brane term on clustering parameter. Furthermore, F(N) is also expressed in terms of b as 
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This is being compared with the earlier distribution function 
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The graphs show the effect of Brane corrections on the distribution function. 



N 


FIG. 1: Comparative study of distribution function F(N) with brane correction and f(N) without brane correction for b = 0.6, 
a = 0.3, iV = 10 and 1 = 1 
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FIG. 2: Comparative study of distribution function F(N) with brane correction and f(N) without brane correction for b = 0.3, 
a = 0.3, iV = 10 and 1 = 1 















V. TWO-POINT FUNCTION 


In this section, we will analyse the power law for galaxy-galaxy correlation function in Brane world. It has been 
demonstrated that the two-point function ^{x) in clustering of galaxies obeys the following power law ® [2^. The 
validity of this power law has been confirmed from N-body simulation [^ . We will now obtain a power law for the 
correlation function using Newtonian potential in brane world models. We will analyse the effect of the the extended 
structure of galaxies and extra dimensions in brane world models on the power law. It will be demonstrated that 
both of these generate higher order corrections in the form of the two-point correlation function. 

In order to obtain the corrections to the power law from the extended structure of galaxies and extra dimensions in 
brane world models, let us consider a spherical system of volume V and energy U, such that it contains N particles. 
Then we can write the following expression [s^, 

U = —NT - - f (j){r)^{r)4Trr^dr. (47) 

2 2 Jy 


Using the interaction potential (j) between two galaxies in brane world models (j) = y 

this expression for the energy of a spherical system (for n = I), we obtain 
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Furthermore, using p = w, we obtain 
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Neglecting higher powers of e/r, we obtain the following power law 
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If we neglect both these effects, we observe that we the power law is very similar Peebles power law ® ~ r~^. 
However, this law receives higher order corrections by taking into account the fact that galaxies have an extended 
structure and by assuming the existence of extra dimensions. The term with e measures the effect of having an 
extended structure, and the term with ki measures the effect of have extra dimensions in brane world models. Both 
these effects have opposite effects on the power law. This it is possible to use this behavior of two-point correlation 
function to test astrophysical data. In case, the astrophysical data shows a positive deviation from the power law, 
it will support the idea that the extended structure of galaxies can not be neglected in analysing the clustering of 
galaxies. However, a negative deviation from the power law will support the existence of extra dimensions in the 
brane world models. 


VI. CONCLUSION 


In this paper, we have analyzed the clustering of galaxies using modified Newtonian potential. The modihcation 
we considered was motivated by the brane world models. It is well known that in brane world models, the standard 
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model fields are fixed on the brane, but the gravitons propagate freely into the bulk. This changes the effective 
gravitational potential, and the deviation from the Newtonian potential can be observed at large distances. Thus, it 
is both interesting and important to use the brane world modified Newtonian potential for analyzing the clustering of 
galaxies. So, in this paper, we analyze the effects of extra dimensions on the clustering of galaxies in the brane world 
models. We calculate the partition functions for a system of galaxies which interact through the modified Newtonian 
potential. This partition function is then used for calculating important the thermodynamic quantities of this system. 
A general clustering parameter is also obtained for galaxies interacting in a brane world model. Finally, we also study 
the effect of the modified Newtonian potential on the two-point function between galaxies. 

We have used the power law deformation of the Newtonian potential which occurs in the Randall-Sundrum model 
where extra dimensions are allowed to be infinite. In other brane world models, the modification of the Newtonian 
potential are parametrized by a Yukawa potential. It would be interesting to repeat this analysis for the Newtonian 
potential modified by a Yukawa potential. It may be noted that the eqnarrays of state for the galaxies interacting 
through the usual Newtonian potential has been obtained and discussed for point mass galaxies [H, [2^. These 
derivations were based on fundamental thermodynamic principles 1^ . Analytic expression for the grand canonical 
partition functions of point masses has also been analyzed [2g. Subsequently, these were generalized to extended 
masses i.e. galaxies with halos, which cluster gravitationally in an expanding universe. From the partition functions, 
the thermodynamic properties, distribution functions, void distribution functions and moments of distributions have 
been obtained more rigorously from statistical mechanical basis [2^. It would be interesting to obtain these results 
for different brane world models and compare the results with observations. It will be interesting to investigate the 
other forms of modified Newton’s law of gravity in other theories like f(R), entropic force gravity and noncommutative 
geometry. We hope to report on these in the future. 
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VII. APPENDIX 


We write the calculations of the first integral of equation (17), the second integral has similar calculations. 
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I^{T,V) = AttV 
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In general 


I,{T,V) = V^ 
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